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ABSTRACT 

Collisions between centimeter- to decimeter-sized dusty bodies are important to understand the 
mechanisms leading to the formation of planetesimals. We thus performed laboratory experiments 
to study the collisional behavior of dust aggregates in this size range at velocities below and around 
the fragmentation threshold. We developed two independent experimental setups with the same goal 
to study the effects of bouncing, fragmentation, and mass transfer in free particle-particle collisions. 
The first setup is an evacuated drop tower with a free-fall height of 1.5 m, providing us with 0.56 s 
of microgravity time so that we observed collisions with velocities between 8 mms -1 and 2 ms -1 . 
The second setup is designed to study the effect of partial fragmentation (when only one of the two 
aggregates is destroyed) and mass transfer in more detail. It allows for the measurement of the 
accretion efficiency as the samples are safely recovered after the encounter. Our results are that for 
very low velocities we found bouncing as could be expected while the fragmentation velocity of 20 
cms -1 was significantly lower than expected. We present the critical energy for disruptive collisions 
<2*, which showed up to be at least two orders of magnitude lower than previous experiments in the 
literature. In the wide range between bouncing and disruptive collisions, only one of the samples 
fragmented in the encounter while the other gained mass. The accretion efficiency in the order of a 
few percent of the particle's mass is depending on the impact velocity and the sample porosity. Our 
results will have consequences for dust evolution models in protoplanetary disks as well as for the 
strength of large, porous planetesimal bodies. 

Subject headings: accretion, accretion disks - methods: laboratory - planets and satellites: formation 
- solar system: formation 



1. INTRODUCTION 



In a protoplanetary disk (PPD) a large number of col- 
lisions between dust agglomerates takes place. The rela- 
tive velocities between these particles, caused by Brown- 
ian motion, differential drift m otions, and gas turbulence 
are depending on their size s ( Weidenschilling & Cuzzi 



1993 Weidling et al.l [2009 , see , e.g. ,~^ According to 



the simulations ol |Zsom et al. (2010) the largest parti- 
cles growing in a minimum mass solar nebula (MMSN) 
model are about two centimeters in diameter. For par- 
ticles of this size, collision velocities of about 10 cms -1 
are expected. However, different PPD models predict 
a wide range of collision velocities for cm-sized dust ag- 
glomerates, r anging from about 1 cms -1 for the model of 
Desch (2007) to about 5 ms" 1 for the low-density model 
By Andrews & Williams (2007). Thus, it is interesting to 
experimentally investigate the collision behavior of cm- 
sized dust aggregates in a wide range of velocitie s also to 
assess the v alidity of the dust-collision model by [Gu ttler 
et al.l ( 2010] ) and the dust-evolution model by Zsom et al"| 
2U1U| ) . For low- velocity collisions between cm-sized dust 
aggregates, it is indispensable to design an experimental 
setup in which the collisions take place without the over- 
whelming influence of gravity. In Sect. [2j we describe 
two experimental setups with which we could investigate 
free collisions between two cm-sized dust aggregates in 
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a wide velocity regime between 8 mms -1 and 2 ms -1 . 
Moreover, the choice of our analog dust particles is ex- 
plained. Sect. [3] shows the results of these experiments 
and Sect. [4]discusses and interprets these results in terms 
of their applicability in PPD dust-evolution models. 

2. EXPERIMENTAL SETUP 

To study the effects of dust aggregate collisions for a 
variety of collision velocities, two slightly different exper- 
iment al se tups were developed: in the drop-tower setup 
(Sect. 2.1 ) we were able to observe free collisions between 
cm-sized dust aggregates for very low to intermediate 
impact velocities (0.8-200 cms -1 ). In a second setup, 
in which we studied velocities above the fragmentation 



threshold (approx. 50 cms -1 , Sect. 2.2), the samples 
were - in contrast to the drop-tower setup - not demol- 
ished after the encounter, so that a detailed measurement 
of the mass balance was possible. We used this setup to 
study the mass transfer between the aggregates as a func- 
tion of the porosity difference and the collision velocity. 

2.1. Drop- Tower Setup 

The first setup is an evacuated drop tower consisting 
of a glass tube with two vacuum chambers attached to 
the top and to the bottom. Outside the glass tube, two 
cameras are recording the experiment in back-light illu- 
mination (Fig. [T]). 

The main part of the drop tower is the glass tube with 
an height of 1.5 m and 22 cm inside diameter. Two vac- 
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Fig. 1. — The drop tower setup consists of (1) a glass tube with 
a height of 1.5 m and an inside diameter of 22 cm, (2) vacuum 
cambers with electrical feedthroughs, pressure gauge and vacuum 
pump, (3) camera platforms with a streamline funnel, (4) two high- 
speed cameras separated by an angle of 90°, (5) buckets filled with 
sand to decelerate the cameras, (6) electromagnets for attaching 
the camera platforms to the ceiling, (7) a two-level particle release 
mechanism (also see Fig. [2j, and (8) two LED arrays and diffusion 
screens for back-light illumination. 

uum chambers are fixed at the top and at the bottom of 
the glass tube, respectively. Electrical feedthroughs are 
attached to these chambers to provide the power for a 
magnetic release mechanism (see below). Additionally, a 
pressure gauge monitors the vacuum quality, which was 
1 — 10 Pa in all experiments. To perform collisions be- 
tween two 2 cm diameter dust aggregates, a two- level 
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Fig. 2. — The two-level release-mechanism: the upper level con- 
sists of a linear solenoid, pressing a nylon string with a metal rod 
against an adjustable metal plate. The lower level is a plate with 
a hole mounted on a rotary solenoid. 



release-mechanism (Fig. [2| is attached at the top flange 
of the upper chamber. The upper level consists of a linear 
solenoid pushing a nylon string embedded in the sample 
onto a vertically mounted plate. The lower level com- 
prises a plate with a hole, mounted to a rotary solenoid, 
where the particle size exceeds the hole diameter and is 
thus centered in position. By moving the position of the 
vertical plate, the impact parameter can be adjusted pre- 
cisely. The distance between the upper and lower particle 
edges can be varied between 1 mm and 30 cm. 

To be able to record a collision in free fall, two high- 
speed cameras are dropped outside the glass tube with 
their fields-of-view centered at the center of mass of the 
free-falling particles. The cameras are each placed on 
a platform mounted on the top of a streamlined funnel 
that ends up with a steel rod. Due to their shape and 
weight of about 10 kg, they are trailing by less than 5 
mm versus the particles (which fall in vacuum) at the 
end of their free fall. Prior to the experiment run, the 
camera platforms are fixed at the ceiling with an electro- 
magnet that can be released by switching off the current. 
The two cameras are separated by an angle of 90° , which 
provides a three dimensional observation of the collision. 
One camera is operated at a picture rate of 261 frames 
per second with a field-of-view of about 6x6 cm 2 and a 
resolution of 190 pixels cm -1 . The other one is a photo 
camera, recording with 40 frames per second and a signif- 
icantly larger field-of-view (15 x 15 cm 2 ) at a comparable 
resolution. A back-light illumination was realized by two 
LED arrays with diffusion screens at the opposite side of 
the glass tube. 

The experiment is triggered by a microcontroller which 
controls two time periods: the first one is the relative de- 
lay between the release of the upper level and the drop 
of the falling cameras. The second one is the time lag 
between the cameras and the lower level. The reachable 
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minimal sum of the two time lags is in the order of one 
millisecond, which leads to the smallest achievable colli- 
sion velocity of 



v co n = g ■ At = 9.81 m s -1 ms ~ 1 cm s 



(i) 



where g is the gravitational acceleration. The upper lim- 
itation of the accomplishable velocities is based on the 
maximal distance between the edges of the particles of 
about Ah = 0.3 m, which yields a collision velocity of 
Vcoii = ^JlgAh w 2.5 ms -1 . 

It is desirable to adjust the collision time to one half of 
the free-fall time, so that the same number of images is 
available for data analysis before and after the collision. 
For this, the distance between the particles Ah must be 
adjusted according to the equation 



Ah v co \\ 



^freefall 



25 



(2) 



with a free fall time for the 1.5 m long drop tower which 

IS tfreefall = 0-55 S. 

The samples used for the particle-particle collisions in 
this experiment are manufactured in a self-built com- 
pression mechanism. It consists of a mold from stainless 
steel and an aluminum stamp with a concave shape at 
one end. The stamp fits into the mold such that a ball 
of 2 cm diameter can be compressed inside. To avoid the 
dust sticking to the mechanism, the surface was worked 
with a nanotechnics polish. For half of the spheres, a 
nylon string is embedded, on which they were suspended 
in the upper level of the particle release mechanism (cf. 
Fig. §. 

For each sample, w e po ured exactly 4.1 g of spherical 
Si02 dust (see Sect. 2.3 for details) into the container 
and compressed it with the stamp. As this pressure could 
not be perfectly identical in all cases, we expect small 
variation in the volume filling factors between <ft = 0.49 
and (j) = 0.54. These values were inferred from the known 
mass and the volume measurement of the falling aggre- 
gates for those cases where the spheres were fully visi- 
ble. Thus, it is not expected that two colliding aggre- 
gates are perfectly similar but one will always be slightly 
more porous and fragile than the other, which will be- 
come important in the discussion. The average difference 
in volume filling factor between two colliding spheres is 
Acj) = 0.025 ± 0.025. 

2.2. Solenoid- Accelerator Setup 

In this series of experiments, collisions between dust 
aggregates with different filling factors at different veloc- 
ities were studied in a velocity range from 0.5 to 1.4 ms -1 
and a range in porosity differences from 0.029 to 0. 129 . 
Dust aggregates of irregular Si02 dust (see Sect. 2.3) 
were prepared in a sim ilar fashion as for the drop-tower 
experiments (Sect. 2.1) by compressing dust into a given 
volume, while we did not aim to produce spherical parti- 
cles here but rather compressed cylinders of 30 mm diam- 
eter and variable heights. We always used the same mass 
of (13.02 ±0.25) g and by qualitatively varying the com- 
pressing pressure we generated dust aggregates, which 
varied in filling factor between <j) = 0.35 and <j) = 0.49. 
This changes the thickness of the cylindrical aggregate 
slightly but typically the cylindrical dust aggregate had 
a height of about 17 mm. In the drop-tower experiments 
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Fig. 3. — The experimental setup for collisions of the dust cylin- 
ders. The solenoid accelerates the lower aggregate, which collides 
with the upper aggregate which is suspended un an elastic strip. 



the aggregates are destroyed after the collision because 
they simply crash on the bottom of the experiment cham- 
ber. Complementary to this, we aimed here to recover 
the aggregates after the experiment to allow mass mea- 
surement, i.e. to determine the mass balance of the ag- 
gregate. To accomplish this, one of the dust aggregates 
was suspended on a thread which itself was fixed to an 
elastic strap. The second agglomerate was launched ver- 
tically upwards and hit the target. The target obtained 
an upward momentum and was somewhat lifted while the 
elastic band decelerated it slow enough on its downward 
motion so that it remained intact and any gained mass 
stayed on the aggregate. With this setup the particle 
could behave in the same way as in a free collision and 
the gained mass stayed intact on the target so this could 
be weighted precisely after the encounter. From first ex- 
periments it turned out that the more compact aggregate 
(higher filling factor) survived the collision and gained 
mass from the less compact aggregate, while this frag- 
mented in the collision. Therefore, the compact aggre- 
gate was suspended and the less compact aggregate was 
accelerated. No alteration of the aggregates' conditions 
could be observed due to the embedding of the string. As 
stable configuration the compact dust aggregate always 
had the same filling factor of 0.480 ± 0.006. The setup is 
sketched in Fig. [3] The launcher consists of a solenoid, 
on which a small platform is mounted. On this platform 
the dust aggregate with higher porosity (i.e. the weaker 
aggregate) is placed. We observed central collisions with 
aggregates colliding with their mantle face. 

The whole assembly was mounted in a vacuum cham- 
ber with an ambient pressure of 0.1 Pa to avoid gas ef- 
fects during the collision. The collisions at the height of 
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the target were captured in an image sequence with 2 
ms time difference between the frames. Based on these 
movies, collision velocities were determined and it was 
assured that the mass gained in a collision remained on 
the target during deceleration. The mass of the target 
was always measured before and after the collisions. 

2.3. Choice of Analog Material 

To prepare samples for an experiment that aims to 
study the formation of planets, proper analog material 
has to be chosen. An additional requirement is the 
commercial availability of the material. We therefore 
used Si02 dust of different shapes in the micrometer-size 
range. This dust has been u sed in many previous experi- 
ments ( |Blum fc Wurm|2008| ) , which showed that it can - 
at least in a mechanical sense - be regarded as represen- 
tative for the class of silicates. For the experiments with 
spherical aggregates in the drop-tower experiments we 
used monomers of 1.5 /am diameter (manufacturer: mi- 
cromod). These g rains have been well characterized by 
Heim et al. ( 1999|) and also the mater ial p roperties were 
analyz ed by Blum & Schrapler (2004) and Guttler et al. 
(2009|. For the cylindrical samples in the second setup 



we used irregular Si02 grains in a size range from 0.1 to 
10 /im where 80 % of the mass are in the range between 1 
and 5 /im (manufacturer: Sigma- Aldrich). Those irregu- 
lar gra ins have also been u sed in many previous experi- 
ments ( |Blum fc Wurm 2008[ ). From these experiments we 



learned that the grain shape plays a minor role for the 
aggregate behavior and those two species can carefully 
be regarded as comparable. We will discuss potential 
differences in Sect. [4] 

3. RESULTS 

In this section, we will present the results of our exper- 
iments, which includes the measurement of coefficients of 
restituti on, s trength of fragmentation and mass transfer. 
In Sect. 3.1 we analyze 24 collisions between spherical 
dust aggregates in a velocity rage of 0.8-37 cms -1 , which 
all led to bouncing. For the higher collision velocities, we 
present eleven collisions of spherical aggregates from 24 



to 202 cms , which are described in Sect. |3.2[ In those 
collisions, one or both aggregates fragmented, and in the 
cases where only one aggregate fragmented, the preserved 
aggregate even gained mass. Above 190 cms -1 , all col- 
lisions led to the disruption of both aggregates. The 
accretion of one aggregate at the intermedi ate collision 
velocity is studied in further detail in Sect. |3.3| There, 



we will present results from collisions between cylindrical 
dust aggregates, in which we varied not only the colli- 
sion velocity but also the porosity difference between the 
aggregates. We performed 24 experiments at a constant 
collision velocity (v = 93 cms -1 ) and varied the porosity 
difference between projectile and target, then 20 experi- 
ments with a constant porosity difference (A</> = 0.091) 
and variable velocities. The particle properties and ex- 
perimental conditions are shown in the table [I] 

3.1. Bouncing Collisions 

In a bouncing collision the two aggregates approach 
each other, lose energy in the dissipative contact, and 
depart with a smaller velocity than before the collision. 
An appropriate parameter to characterize these collisions 




Fig. 4. — A typical sequence of a bouncing collision at an impact 
velocity of v = 8 mms - 1 . The coefficient of restitution for this 
experiment was e = 0.29. 



is the coefficient of restitution 

Uafte 



e = 



^befo 



(3) 



The relation 1 — s 2 is a measure for the energy dissipa- 
tion. Figure [4] shows a typical sequence of a bouncing 
collision, where the impact velocity was v = 8 mms -1 . 
The second image is close to the moment of contact and 
already these raw images at constant time interval allow 
a good guess for the coefficient of restitution (the exact 
value was calculated to e = 0.29). 

The main parameters to measure in these collisions are 
the velocities before and after the collision. To achieve 
this, we binarized the image sequences of both cameras 
and convolved them with a disk of the proper sphere 
diameter to compute the temporally resolved, three di- 
mensional position of the two aggregates. This yielded 
a precision of 2 pixels, therefore the quality of the three 
dimensional distance between the center of the two parti- 
cles is less than 0.2 mm for a single frame, which becomes 
better when using a sequence of images. The distance 
was separately fitted in all three dimensions (for the ver- 
tical direction we chose the fast camera) to get a velocity 
vector before and after the collision. From these, we can 
calculated the coefficient of restitution e, the normalized 
impact parameter b/R, and the absolute impact velocity 
v. A rotation of the lower particle was not found in any 
case. In some cases the upper particle suspended at the 
embedded nylon string librated, which led to a rotation of 
the particle in the free fall. To measure the rotation en- 
ergy and the surface velocity, the position of two points 
on the surface relative to the center were followed and 
the surface velocity and the rotation energy calculated. 
That energy was compared with the translational energy 
and a threshold was set at 10 percent to distinguish be- 
tween rotating and non rotating impacts. Below, we will 
consider only central collisions with impact parameters 
b/R < 0.3, and experiments with rotating samples are 
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TABLE 1 

Overview on the conducted experiments and their outcomes 
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Fig. 5. — The coefficient of restitution e does neither show a 
correlation with the impact parameter (top) nor with the impact 
velocity (bottom). A correlation analysis yields values of —0.003 
and 0.1, respectively. Triangles mark all three dimensional data 
and additionally the 2 two dimensional data are shown (bottom) 
as squares. 

neglected as well. 

Let us first evaluate the dependance of the coefficient 
of restitution on the normalized impact parameter b/R. 
As shown in Fig. [5] (top), no strong correlation between 
these two parameters exists. However, to verify the in- 
dependence, an analysis for the correlation of these pa- 
rameters was performed. With this correlation analysis 
the parameters b/R and e could clearly be determined 
as uncorrelated in the considered parameter range. In 
Fig. [5] (bottom) the coefficient of restitution is plotted 
over the impact velocity. Here, we increased the number 
of data points by also including 7 collisions where only 
the high-speed camera was operative, thus where we only 
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Fig. 6. — The normalized cumulative number over the coeffi- 
cient of restitution e suggests an integrated Gaussian function what 
yields to a normal distribution for e. The average is at 0.35 ±0.12. 
The dashed and the dash-dotted lines mark one and two standard 
deviations, respectively. 

have two dimensional information. The collision veloc- 
ity is dominated by the known vertical component and 
the coefficient of restitution was corrected by using the 
average of the horizontal dimension of the photo-camera 
from the three dimensional cases. This led to a factor 
of 0.036 that was added to the coefficient of restitution 
for these collisions. The error bars are as big as the 
range of the photo-camera component finding in the 3 
dimensional collisions. With these 24 collisions, Fig. [5] 
(bottom) does not show an obvious correlation between 
impact velocity and coefficient of restitution and to verify 
this, we again performed a correlation analysis of these 
parameters which yielded a correlation coefficient of 0.1. 

Conclusively, all impacts can be regarded as cental or 
near-central collisions and the coefficient of restitution is 
neither correlated with the impact parameter nor with 
the collision velocity in the parameter range considered. 
We thus plotted the cumulative number fraction over the 
coefficient of restitution (Fig. [6|. The fit is an integrated 
gaussian function 



P(e) 



f — 
Jo cr £ V27r 



exp 



de' 



(4) 



with an average coefficient of restitution of e = 0.35 and 
a standard deviation of a £ = 0.12. The gray dashed 
and the dash-dotted lines mark one and two standard 
deviations, respectively. 

3.2. Fragmenting Collisions 

Increasing the impact velocity led to aggregate frag- 
mentation. In these cases, the coefficient of restitution 
could not be calculated, but the strength of fragmenta- 
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Fig. 7. — An image sequence of a fragmenting collision at a ve- 
locity of 47 cms -1 . The upper aggregate is destroyed while the 
lower one gains 0.6 percent of its original volume. 

tion fi is the characterizing parameter which we deter- 
mined. It is defined as /ij = M^^/Mq^ were Mf^ is the 
largest fragment and Mq^ is the total mass of the aggre- 
gate i (with i = 1 or 2) before fragmentation. The colli- 
sion velocity was determined in the same way as in Sect. 
|3.1| and the size of the largest fragment was measured 
on one or more representative images on the photo cam- 
era. Again, these images were binarized and the cross- 
sectional areas of the fragments were fitted with ellipses 
of the same expanse. A lower estimate for the fragment 
volume follows from an ellipsoid where the ellipse is ro- 
tated around the long axis and the upper estimate is 
for rotation around the short axis. Our best estimate is 
the mean value between those two and the error follows 
from the values themselves. To calculate the mass of 
the fragments we assumed that the volume filling factor 
remained unchanged. The results on the fragmentation 
strength (i.e. the relation fifv)) will be discussed in Sect. 

i 

The slowest impact velocity where fragmentation oc- 
curred was at 24 cm s _1 and the fastest bouncing collision 
was at 37 cms -1 . In Fig. j\ a sequence of a collision at 
47 cms -1 is presented. This is a typical collision if the 
velocity is below 187 cms -1 , as for these velocities only 
one aggregate fragmented while the other one was intact 
and even gained mass. This roughly conically shaped 
dust pile is marked with a green circle in the last image 
of Fig. [7[ We assume this pile is rotational symmetric, 
we could easily compute its volume: we rotated the im- 
age such that the long axis was vertical and measured 
the width of the pile for each line. Due to the symme- 
try, each line represented a disk with a height of one 
pixel and by summing up these disks we got a volume. 
Again, we assumed that the volume filling factor was 
unchanged to get the mass of the dust pile. In the ex- 
ample in Fig. [7] the intact aggregate has grown by 0.57 
percent of its own mass. The mass gain for all 6 exper- 
iments ranges from 0.05 to 1.1 percent and the velocity 
dependence will be presented in Sect. |3.3| There was 
no relation on whether the upper or lower particle frag- 
mented as these two possibilities exactly split up to the 
6 experiments. However, the effect that only one parti- 
cle was fragmented had a clear limitation at a maximum 
velocity of about 190 cms -1 . The only exception was 
one collision at v = 47 cms -1 , in which both aggregates 
fragmented. For all collisions above 190 cms -1 both par- 
ticles fragmented. 

Depending on collision velocity, we qualitatively found 
the three different collisional outcomes that comprise 
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Fig. 8. — The strength of fragmentation /i over the impact veloc- 
ity for all drop-tower collisions. The axis for fi > 1 was stretched for 
better visibility, /j is plotted for each individual sphere so that we 
get two data points for each experiment. The values for fragmenta- 
tion (red) are fitted with a power law (dash-dotted line) and those 
for mass gain (green) linearly. The black dashed line at fi = 0.5 
denotes the limit for catastrophic fragmentation (cf. Sect. El. 

bouncing, fragmentation with mass transfer and com- 
plete fragmentation. In Fig. [8] we present a conclusive 
overview of the outcomes in the drop-tower experiments. 
The strength of fragmentation for each individual sphere 
in a collision fii is plotted as a function of the collision 
velocity. For the bouncing collision (previous section) 
the mass remains constant, so that the value for these 
collisions is unity. Each fragmenting collision yields two 
H values, i.e. one for each aggregate. The disruptive 
collisions are given by the red squares and those which 
lead to mass growth yield one negative and one positive 
/i value (red and green triangles) for the disrupted and 
grown aggregate, respectively. We adopt a power law to 
the values with \i < 1, which is given by 



-(") 



(0.18 ±0.04) ms- 



-0.52±0.09 



(5) 



and shown as a red dash-dotted line. The vertical axis 
for the collisions with mass gain is linear and stretched 
and we found that it can be described by a linear relation 

ii+{v) = (-4.3 • nr 4 + 6.o • nr 3 • ( — ± 1. (6) 

Vms 1 / 

More details on the accretion efficiency will be given in 
the next section. 

3.3. Accretion Efficiency 

To further study the accretion efficiency as a function 
of the porosity difference and collision velocity, we first 
kept the porosity difference constant at 0.091 ±0.014 for 
20 experiments and varied the collision velocity from 0.5 
to 1.4 ms _1 . We used the setup shown in Fig. [5] for 
this. In a second step, we performed 24 experiments at 
a constant collision velocity of (0.93 ± 0.08) ms" 1 and 
varied the porosity difference between the aggregates in 
the range of 0.03 to 0.13. The target aggregate - the 
upper aggregate, which was more compact and thus ro- 
bust - always gained mass as can be seen in Fig. [9] and 
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t = 1 74 ms 



Fig. 9. — An image sequence that illustrates a collision between 
the cylindrical samples. Here, the collision velocity is roughly 
1 m s ~ 1 . 

after capturing the intact aggregate, we were able to mea- 
sure the accretion efficiency e ac , i.e. the mass growth in 
units of the projectile mass. For a few examples, the 
mass growth was also deduced from the images with the 
same approach as presented in the drop-tower experi- 
ments. Here the deduction of mass growth from the im- 
ages are in good agreement with the measured accretion 
efficiency. For this reason it can be concluded that the 
volume filling factor of the cone remains unchanged in 
comparison with the projectile filling factor. 
The mass gain dependence on the collision velocity is 



presented in Fig. 10 , the gray circles represent single 
experiments. Below 0.5 ms -1 no mass growth but only 
bouncing collisions without mass transfer were observed 
(two data points at e ac = 0). Therefore, we can confirm 
a low-velocity cut off, where collisions are not energetic 
enough for sufficient inelastic interaction, i.e. to lead to 
fragmentation. The values above 0.5 ms -1 are averaged 
in steps of 0.2 ms -1 up to 1.5 ms -1 , which is shown 
by the black squares. An increase of the mass gain with 
collision velocity is clearly visible. Faster collisions lead 
to a higher mass gain. We assumed a linear increase and 
fitted a straight line 

e ac (v, = 0.091) = -3.5 • 10~ 3 + 2.5 • 10~ 2 ( —^A 

Vms _i / 

(7) 

to the raw data (neglecting the first point at e ac = 0). 
The black triangles are the results from the experiments 
in the drop-tower setup presented in the previous section 
and show a similar trend. The lower values can be ex- 
plained by a smaller difference in porosity (see below and 
Sect. [§. 
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Fig. 10. — Mass gain over collision velocity at constant porosity 
difference Acf) = 0.091. Each gray circle represents one experiment, 
the black line a linear fit to these data points and black squares 
are averaged values to illustrate the linear trend. Results from 
the drop-tower are presented as triangles and are also linearly fit- 
ted (lower solid line). The dotted lines are the result from the 2 
dimensional fit which is discussed in Sect. @] 
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Fig. 11. — Mass gain as a function of porosity difference for a 
constant collision velocity of 0.93 ms -1 . The gray circles show the 
raw data, the solid curve is a fit to these data according to Eqs. [8] 
and[9] these equations are undistinguishable from each other in the 
parameter range shown, and the black squares represent averaged 
values. The triangle shows the mean accretion efficiency of all drop 
tower experiments. The dashed curve is discussed in Sect. ^] 



The mass gain as a function of the measured poros- 
ity difference of colliding aggregates is shown in Fig. [TTJ 
again, the gray circles denote individual experim.ents.~~A 
trend of increasing accretion efficiency with increasing 
porosity difference is already visible in the raw data and 
becomes even clearer after averaging in steps of 0.02 in 
porosity. Without further justification but to give a sim- 
ple mathematical expression we fitted a parabolic curve 
to the raw data given as 

e a c,par (A0, v = 0.93 ms -1 ) = 3.4-lO" 3 +1.44-A0 2 , (8) 
which is shown as a solid line in Fig. [TT] Alternatively, 
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Fig. 12. — The coefficient of restitution of various experiments 
as a function of the volume filling factor, although very clear, 
this trend should only be regarded as qualitative because of the 
difference in samples, impact parameter and collision velocity. 

an exponential function given as 



e acexp (A</),v = 0.93 ms -1 ) = 4.4-10~ 3 +1.2-10~ 3 -exp(23 

(9) 

can hardly be distinguished from the previous curve and 
gives a similar (even slightly better) reduced x 2 . Both 
curves fit the data in the presented range while they 
significantly deviate for higher porosity differences: the 
exponential function gives an accretion efficiency of 0.5 
for a porosity difference of 0.26, while for the parabolic 
function this happens at a porosity difference of 0.59. 
An extrapolation of this curve is dangerous and not 
recommended but the consequences would be that the 
parabolic function predicts a limeted mass transfer even 
for the highest possible porosity differences while the ex- 
ponential function predicts a saturation or the onset of 
a new effect for porosity differences as small as 0.3. This 
can be tested in future experiments. 

4. DISCUSSION 

Depending on collision velocity, we qualitatively found 
three different collisional outcomes that comprise bounc- 
ing, fragmentation with mass transfer and complete 
fragmentation. In Fig. [8] we presented a conclusive 
overview of the outcomes in the drop-tower experiments, 
which give us the velocity range of these effects: 



Fig. 13. — Two-dimensional fit interpolating the mass accretion 
as a f unc tion of velocity and porosity difference according to Equa- 
tion ( |10| . The data points represent the measured accretion effi- 
ciency m color, for the experiment with cylinders we only show the 
mean values. 

ments and a similar velocity (20 cms -1 ), the aggregates 
A^)had a cubic shape and the volume filling factor was much 
lower ((/> = 0.15). |Blum fc Miinch| ( |1993[ ) performed near 
central collisions between aggregates of ZrSi04 dust with 
a comparable grain size and a volume filli ng factor of 
= 1.26. The fourth data point is from HeiBelmann 



v < 40 cms 1 
v > 20 cms -1 and 
v < 190 cms -1 
v > 190 cms -1 



et al. (2010) who used solid glass beads (still Si02) with 
all possible impact parameters in a multiple collision ex- 
periment within a range of velocities (0.3 to 5 cms -1 ). 
In spite of these differences, a clear trend is obvious, 
which is that porous aggregates behave more inelastic, 
thus dissipate more energy than compact aggregates or 
even solid samples. This trend is much stronger than the 
dependence on the collision velocity, which we could not 
resolve in Sect. 13.11 

For velocities around 20 cms -1 we find that particle 
fragmentation sets in, which means that only one aggre- 
gate is destroyed and transfers mass to the other par- 
ticle. In Figs. [lO] and |TT] we presented the accretion 
efficiency as a function ofcollision velocity and porosity 
difference. From these, we now want to give a full equa- 
tion e ac (v,A(j)). We cannot simply multiply Equations 
Q and (or (J9|, respectively) because these do not 
exactly match at v = 0.93 ms -1 and A<j) = 0.091. We 
therefore fitted the raw data points from the experiments 
with the cylindrical samples with the function 



bouncing 

partly fragmentation with 
mass transfer 
disruptive fragmentation 



e ac (v, A<j>) = (a + b • • (c + d • v) 



(10) 



Bouncing collisions were presented in Sect. 3.1 and 
the main result was the mean coefficient of restitution 
of £ = 0.35 ±0.12, where the error does not represent 
the measurement precision but the standard deviation 
resulting from the natural scatter of the data. To com- 
pare this results with earlier experiments, we present the 
coefficient of restitution as a function of the volume filling 
factor for various experiments (Fig. [l2|). These exper- 



A special search technique like the downhill simplex al- 
gorithm was used to obtain the parameters. The calcu- 
lated parameters are a = 0.26, b = 85, c = — 6.1-10 -3 and 
d = 2.5 • 10 -2 s m -1 . Equation ( 10 ) can again be reduced 



iments are all different: while HeiBelmann et al. (2007) 
used the same dust material as in our drop-tower experi- 



to the constant porosity difference and constant veloc 
ity for which the series of experiments were conducted, 
which is presented in Figs. [lOland 11 as the dashed line. 
Those are still in the range of the errors and differ only 
slightly from the original functions. The combination of 
both parameters is presented in the contour plot of Fig. 
13 The color values correspond to Equation (TTOJ) and 
the color of the square symbols represents the value of 
the averaged data (squares in Figs. 10 and 11). The 
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Fig. 14. — Fragmentation strength as a function of impact ve- 
locity. The triangles, squares, a nd ci rcles show the d ata from Fig. 
M from |Blum & Munchj (1993), and |LammeT[ ( |2008| >, respectively. 
Th e dashed and solid lines are power-law hts to these data (see Eq. 
[Rf. 



white region at low velocities (v < 25 cms -1 ) is where 
Equation (10) mathematically predicts a negative mass 
transfer, while we rather expect bouncing collisions. This 
is supported by the black triangles, representing bounc- 
ing collisions without mass transfer. Due to the difference 
in aggregate geometry and grain shape, the results from 
the drop-tower experiments were not utilized to compute 
the fit in Equation ( flO] ) but they are nonetheless very 
well represented by this curve (see triangles in Fig. 13). 
If we again reduce it to the constant porosity difference 
= 0.025, we find in Fig. [To] that the dashed line 
is very close to the solid line. This agreement without 
further assumption assures us that the experiments are 
highly comparable in spite of the discussed differences. 
We therefore regard Equation ( 10 ) as a good description 
for the mass transfer in a disruptive collision between two 
aggregates with different strength. The validity range of 



this equation is given by 20 cms - < v < 190 cms 
in velocity. For the porosity, we cannot make a definite 
prediction for A<p > 0.13, which is the maximum value 
that was studied here. 

The velocity of 20 cms -1 where the fragmentation sets 
in is significantly lower than expe cted from experiment s 
with millimeter-sized aggregates (Guttler et al. 2010). 
We therefore expect a size dependence ol the t hresh- 
old velocity and com pare our results with data of [Blum 
|fc Mimch| (jl993[ ) and |Lammel| ( |2008D who measured the 
strength ol fragment ation lor a variable velocity as pre- 
sented here in Sect. 13.21 We consider central collisions 



and ta ke the mean values from Fig. 7 i n Blum fc Munch 
( |1993| ) as well as raw data of |Lammel| ( |2008p . This is 
presented in Fig [14] w here the red triangles represent our 
da ta from Sect. |3.2| the green squares are the values 
of Blum & Miincn (1993, mean values, thus including 
erro r s) and the Blue circles represent the data of Lam- 



mel (2008). From this we get equations similar to Eq. 
([5j, which relate the fragmentation strength to the colli- 
sion velocity between two dust aggregates. In its general 



Fig. 15. — The critical fragmentation strength Q* as a function of 
dust-aggregate radius for o ur data (triangles) as we l l as fo r the data 
of |Blum fc Munch] (1993) (squares) a nd |Lammel| ( |2008]> (circles). 
A least-squares power law fit after Eq. |14| is shown as solid line for 
the two cases u = 1.0 a nd u = 0.5. Also shown are the Q* data of 

(crosses) as well as the 
line). The dotted line 




Asphaug ( 1999 ) 



form, this reads 



according to the model by |Benz 



(ii) 



Our data, presented in Fig. 



yield v 



(Uif-o.oD ms- 1 and a = 0.52 ± 0.10 for collisions be- 
tween equal-sized dust aggregates with masses of m = 
4.1 g and radii of s = 1 cm (short dashed red line in 
Fig. " 

gates 



14) 



For collisions between mm-sized dust aggre- 
l.e. for the data of |Blum fc Munch] ( | 19931 ) (using 
aggregates with m = 4 mg and s = 0.9 mm consisting of 
/xm-sized ZrSiC^), we get = (1.01 ±0.03) ms -1 and 
a = 0.38 ± 0.04 ( long dashed, green line). For the data 
of Lammel ( 2008 ) (using aggregates with m = 1 mg and 
s = 0.7 mm consisting of /im-sized Si02; s is the half edge 
length of the cubic aggregates here), the parameters are 
v = (1.73±g;||) ms" 1 and a = 0.54 ± 0.11 (solid blue 
line). These relations between fragmentation strength 
and velocity can be used to derive the velocity for the 
onset of fragmentation, i.e. i>i.o = v(fi = 1.0) = vo 
and for the case that the largest fragment possesses 
half the mass of the original projectile aggregate, i.e. 
^0.5 = v([i = 0.5) = vo • 0.5 -1 /". Applying this to all 
the data in Fig. [TZJ we get ^0.5 = 0. 64 ms -1 for our 
data, ^o.5 = 6.2 ms -1 for the data of |Blum & Munch 



(1993) as well as for the data of Lammel|" ( 2008| ); the ve- 
locity for the onset of fragmentation is already given by 
the fit parameters vq. These velocities can be converted 
into specific energies, assuming equal-mass collision part- 
ners, by the relation Q x = ^v%, with x = 1.0 or x = 0.5. 
We want to bring forward a relation for the critical frag- 
mentation energy depending on aggregate size but have 
to be careful in comparing these data as the volume fill- 
ing factors significantly differ from 0.15 (Lammel 2008) 
to 0.5 (this work). 

The porosity has an effect on the strength of the ag- 
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gregates as well as on the transmission of impact stress 
waves. The tensile strength for these aggregates is given 

by 

T(0) = -(lo 2 -^ 1 - 48 ^) (12) 



(Guttler et al. 2009| their Eq. 8) and for the relation 



between aggregate strength and fragmentation energy we 
consider the r esults of |Setoh et al. ( |2010| ) and |Machii| 
' fc Nakamura] ( |201l[ ) and take Q oc T 3/4 . The effect of 
poor shock transmission in porous media was studied by 
Love et aT] ([1993) who found a porosity dependence of 
Q oc (j)~ 6 '° if the strength of the bodies remains constant 
for changing porosity. In their numer ical simulation on 
the disruption of larger bodies, also Jutzi et al. (2010) 
found that por ous targets are m ore difficult to disrupt. 
Thus, following Love et al. (1993), we apply a correction 



Q*(0 = O.5) = Q a 



3. 61 nl- H-0.5 



0.5- Jb 10 



-3.6 10 1 



(13) 



where Q x {<j>) is the critical energy for the individual ex- 
periments as calculated above and Q x (4> = 0.5) is the 
critical energy corrected to a volume filling factor of 
(j) = 0.5 as used in our experiments. With this, we ar- 
rive at Fig. [l5j where the critical fragmentation energies 
Qo.5 and Qi.o for the three experiments are presented 
as a function of aggregate size. A power-law fit to these 
data of the form 



-P 



J kg" 1 



(14) 



yields f3 = 0.95 ± 0.38, 
si.o = (3.4t?1) • 10- 5 



-1.2) ' 10 " 5 m - We have to note that the 
slope should be regarded as an estimate and it is de- 
sirable to have another experiment to measure the frag- 
mentation energy for different sized samples but with the 
same strength and porosity. Accordant studies are cur- 
rently on the way. Our slope is slightly ste eper than in 
the material- strength regime presented by |Benz fc As 



phaug (1999) (see Fig. 15), who studied collisions be- 



tween rocky bodies in computer simulations. Our firm 
absolute values for Q0.5 are several orders of mag nitude 
smaller than tho se experimentally measured by J Set oh] 
et al.| ( |2007D and |Ryan et~aT1 \ 1991 1 ) (see Fig. [15} , who 
used sinter ed glass beads and w e akly p lued gravel, re- 
spectively. Machii & Nakamura (2011 
experiments with sintered glass beads (I 
Fig. 15) and found a dynamic strength in t he same or 
der oimagnitude as Setoh et al. (|2007|) and Ryan et al. 



also performed 
lack diamond in 



(1991). However, they also measured the static tensile 
strength of the material (3.5 • 10 5 Pa) and if we apply 
a similar scaling as above, we arrive at the gray dia- 
mond which is in good agreement with our results. Al- 



though not p erfectly applicable f or our samples, we use 
the model by Benz & Asphaug (1999) for the strength 
in the gravitational regime. Actually, in a recent study 
of Jutzi et aT] (|2010|) the slopes in the gravity regi me (cf. 



their Tab l e 3) a re mostly comparable to those of 



Benz 



Asphaug] ( |1999| ), also if they use porous bodies, so that 
we take 



^0.5 



1 mm 



-0.95 



(15) 



illustrated by the dotted line in Fig. [15] The first 
term is simply given by Eq. [l4| while the second term is 



taken from the gravity regime in |Benz fc Asphaug (1999 
cf. their Eq. 6 and Table 3) assuming a bulk density 
of 1000 kg m -3 . Thus, our measurements suggest that 
the weakest bodies in PPDs have radii of the order of 30 
cm. Their strength should be around Q* = 10 -2 J kg -1 , 
which corresponds to impact velocities between equal- 
sized particles of about 0.2 ms _1 . A comparison to the 
expected collisi on velocities in differe nt PPD models 
(see Fig. 8 in Weidling et al. 2009[ ) shows that the 
collision velocities lor decimeter-sized aggregates ex- 
ceed this critical fragmentation veloci ty in both, the 
minimum-mass solar nebula model ( Weidenschilling 



1977[) as wel l as in the low-density model by | Andrews 
fc Williams| (2007). The collision velocities in the 



high- density model by |Desch| ( |2007[ ) are comparable to 
our estimated fragmentation velocity for decimeter-sized 
dust aggregates. Thus, dust aggregates of decimeter 
size can only barely survive mutual collisions if the 
impact velocities are very moderate. Alternatively, the 
strength of aggregates against fragmentation increases if 
the aggregates experience any solidification p roces s, e.g. 
by sintering (see data by [Setoh et al. ( |2007[ ) and [Ryan 
[etlnT] ( [l99lf in Fig. [T^). 
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